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Abstract 

We discuss how to reduce theoretical uncertainties in the neutrino-deuteron breakup cross- 
sections crucial to the Sudbury Neutrino Observatory's efforts to measure the solar neutrino flux. 
In effective field theory, the dominant uncertainties in all neutrino-deuteron reactions can be ex- 
pressed through a single, common, isovector axial two-body current parameterized by L\^a- After 
J> . briefly reviewing the status of fixing L% a experimentally, we present a constraint on L\ a imposed 

by existing reactor antineutrino-deuteron breakup data. This constraint alone leads to an uncer- 
tainty of 6-7% at 7 MeV neutrino energy in the cross-sections relevant to the Sudbury Neutrino 
Observatory. However, more significantly for the Sudbury experiment, the constraint implies an 
uncertainty of only 0.7% in the ratio of charged to neutral current cross-sections used to verify 
the existence of neutrino oscillations, at the same energy. This is the only direct experimental 
constraint from the two-body system, to date, of the uncertainty in these cross-sections. 
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I. INTRODUCTION 



Recent results from the Sudbury Neutrino Observatory (SNO) [|TJ highlight the impor- 
tance of a precise determination of neutrino-deuteron breakup reaction cross-sections. The 
three reactions used by SNO to detect the 8 B solar flux are 

u e + d p + p + e~ (CC), 
v x + d — > p + n + v x (NC) , 
v x + e"^ v x + e~ (ES). 

The charged current reaction (CC) is sensitive exclusively to electron-type neutrinos, while 
the neutral current reaction (NC) is equally sensitive to all active neutrino flavors (x = 
e,n,r). The elastic scattering reaction (ES) is sensitive to all active flavors as well, but 
with reduced sensitivity to z/ M and v T . Detection of these three reactions allows SNO to 
determine the electron and non-electron active neutrino components of the solar flux, and 
it is then obvious that the cross sections for these three reactions are important inputs for 
SNO. However, while the ES cross section is very well determined, the CC and NC cross 
sections have never been tested to high precision. 

Theoretically, the complications in the CC and NC processes are due to two-body currents 
which are irreducible interactions involving leptonic external currents and two nucleons. 
The two-body currents can be calculated either through meson exchange diagrams aided 
by modeling of any unknown weak couplings, or can be parameterized using effective field 
theory (EFT). In both cases, experimental data from some other process are required in 
order to calibrate the unknowns in the problem. In EFT, this calibration procedure can be 
described in an economic and systematic way. The reason is that, up to next-to-next-to- 
leading order (NNLO) in EFT, all low-energy weak interaction deuteron breakup processes 
depend on a common isovector axial two-body current, parameterized by L^a 0- This 
implies that a measurement of any one of the breakup processes could be used to fix L^a- 

In this paper, we will briefly review the EFT approach and discuss experiments that could 
be used to fix L\,a- Then we present the constraint on L^a using reactor Vd scattering. 



II. EFFECTIVE FIELD THEORY 

For the deuteron breakup processes used to detect solar neutrinos, where E u < 15 MeV, 
the typical momentum scales in the problem are much smaller than the pion mass m n (~ 140 
MeV). In these systems pions do not need to be treated as dynamical particles since they only 
propagate over distances ~ l/m,, much shorter than the scale set by the typical momentum 
of the problem. Thus the pionless nuclear effective field theory, EFT(f) || |], [5|, ||, 0, is 
applicable. 

In EFT (ft), the dynamical degrees of freedom are nucleons and non-hadronic external 
currents. Massive hadronic excitations such as pions and the delta resonance are "integrated 
out," resulting in contact interactions between nucleons. Nucleon-nucleon interactions are 
calculated perturbatively with the small expansion parameter 

Q= (1/Q A 7 ' P) (1) 

which is the ratio of the light to heavy scales. The light scales include the inverse S-wave 
nucleon-nucleon scattering length l/a(^ 12 MeV) in the 1 5'o channel, the deuteron binding 
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momentum 7(= 45.7 MeV) in the 3 Si channel, and the typical nucleon momentum p in 
the center-of-mass frame. The heavy scale A is set by the pion mass m n . This formalism 
has been applied successfully to many processes involving the deuteron || |9|, including 
Compton scattering |L0], [□]], np — > for big-bang nucleosynthesis [12|, jl3"| , vd scattering 
for SNO physics 0], the solar pp fusion process \TA, and parity violating observables |TB 



Also studies on three-nucleon systems [I7| have revealed highly non-trivial renormalizations 
associated with three body forces in the S\/2 channel (e.g., 3 He and the triton). For other 
channels, precision calculations were carried out to higher orders |7], 

In Ref. ||, EFT(7f) is applied to compute the cross-sections for four channels (CC, NC, 
V e + d^e + + n + n (77CC) and V x + d — > V x + n+p (z7NC)) to next-to-next-to- leading order 
(NNLO), up to 20 MeV (anti)neutrino energies. As already mentioned, these processes have 
been shown to depend on only one parameter, L\ a, an isovector axial two-body current. 
This dependence is subject to an intrinsic uncertainty in our EFT calculation at NNLO of 
less than 3%. Through varying L^a, the potential model results of Refs. 
reproduced to high accuracy for all four channels. This confirms that the 
between Refs. 



are 



1J and |20 
~ 5% difference 



T9| and is due largely to different assumptions made about short distance 



physics. 

The same two-body current L\^a also contributes to the proton-proton fusion process 
p + p ^ d + e + + u e . This is the primary reaction in the pp chain of nuclear reactions that 
power the sun, reactions which in turn generate the neutrino flux to be observed by SNO. 
The calculations in EFT(/f) were carried out first to second order 
order 1X31 . 



131, and then to fifth 



III. FIXING L 1>A 

In order to determine neutrino-deuteron breakup reaction cross-sections and the pp fusion 
amplitude to high precision, one needs a precise determination of L^a- Naive dimensional 
analysis gives |Li,a| ~ 6 fm 3 when the renormalization scale \x is set to m^. This implies that 
the contribution from L^a is & t the 7-8% level at a neutrino energy of 7 MeV. However, 
this is a /i-dependent statement and cannot be used to draw comparisons between two- 
body contributions in the EFT calculation and conventional potential model calculations. 
It is therefore important to consider how one can constrain L\ t A experimentally. This is a 
daunting task, as most weak processes where L^a contributes are difficult (if not practically 
impossible) to measure accurately in the laboratory. Further, some experiments such as 
the measurement of the flux-averaged CC cross-section using neutrinos from stopped muon 
decays J2T], employ neutrino energies greater than 20 MeV. At this time, the convergence of 
the calculations of Ref. is uncertain at these higher energies, so we are forced to ignore 
such experiments. Here we list some observables that could be used to determine L\ a- 



1. v e + d — > e~ + p + p: the planned ORLaND detector |22| has proposed to measure 
this CC process with ~ 1% accuracy. This measurement, combined with higher order 
calculations in EFT, would calibrate SNO's CC and NC processes to the same level 
of accuracy. 

2. reactor V e d breakup reactions: the main topic of this letter. We discuss the extraction 
of L\ a in next section. 
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3. solar v x e elastic scattering (ES): the three channels available for measuring the solar 
neutrino flux, ES, CC, and NC are all measured by SNO, and ES is measured at Super- 
PC to high precision. In general, these results can be used to constrain three quantities: 
the electron and non-electron active neutrino components of the solar flux, and L\^. 
As statistics improve in all three channels, SNO could become self-calibrating in that 
it could use this determination of L\ a to fix the CC and NC cross-sections |213 . 



4. muonic atom capture /i~ + d — > u^ + n + n |2"4] , ^J: this reaction can involve significant 
energy and momentum transfers. EFT (7^) may fail in regions of phase space where the 
neutrons are energetic, but should be rapidly convergent in the region of phase space 
where the two neutrons move slowly [p5| . Coincident measurements of those slow 



neutrons to the required precision are very difficult. However, indirect high-precision 
measurements of the total capture rate are feasible at PSI by comparing the lifetime of 
\i~ to that of /i + in a deuterium target (thereby avoiding the need to detect final-state 
neutrons), then subtracting the (easier to measure) faster neutron events p6 . 



5. Tritium beta decay H — > He + e + + v e : under the assumption that three-body 



currents are negligible, Schiavilla et al. |27j used this process to fix the two-body 



current and made a prediction for the solar fusion process p + p — > d + e + + v e , with an 
accuracy of better than 1%. This prediction by Schiavilla et al. [|7fl can be translated 
to a constraint on L^a using the EFT formula of [TJ]. After updating the pion-nucleon 



coupling g A from 1.26 to 1.267, we obtain Lg iaviUa et aL = 6.5 ± 2.4 fm 3 at NNNLO. 
This result includes the theoretical uncertainty discussed in Ref. l~5j| . If we truncate 



the results from Ref. [D| at NNLO, the expression for A(0) becomes 



Annlo(O) = 2.61 + 0.0104^-^4 



^tvlo(O) = 2.61 + 0.0104(— + 0(1.5%) (2) 

If we compare this to Schiavialla et al., we find ^Schiavilla et ai. _ 4 2 ± 3 7 f m 3 ^his 
latter result is more appropriate for comparison to the constraints presented later 
in this work. We note also that, following the approach of Schiavilla et al., other 
calculations have made use of tritium beta decay to make predictions for other weak 
processes, and all claim an accuracy of ~ 1% 



6. Helioseismology: the ability of the standard solar model to reproduce acoustic mode 
(p-mode) oscillations in the Sun to high precision can be used to constrain the pp 
fusion cross-section, and thus L^a- Assuming that the solar model must reproduce 
these oscillations to the same precision, it is found that L M = 7.0 ± 5.9 fm 3 gjj at 



NNNLO, or Li t A = 4.8 ± 6.7 fm 3 using eq. ^| at NNLO. These numbers also include 
the theoretical uncertainty in the reduced matrix element for pp fusion |15|] (but not 
from other theoretical uncertainties in the solar model). It should be noted that this 
constraint is more an indicator of the scale and sign of There are other physical 
inputs to the solar model (e.g., opacities) that will weaken this constraint if also allowed 
to vary. 

So far, we have only given numerical constraints on L\ a i n methods 5 and 6. These 
constraints rely on certain assumptions that require deeper theoretical study. In the next 
section we will look at fixing L 1; a using method 2, namely reactor V e d breakup reactions. 
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IV. A CONSTRAINT FROM REACTOR ANTINEUTRINO EXPERIMENTS 



The charged current V e d — » e + nn (fCC) and neutral current V x d — > V x np (z/NC) deuteron 
disintegration have been observed in several experiments with reactor antineutrinos. We can 
use the results of these experiments to constrain the parameter Lxa- 



The results of the pioneering experiment at Savannah River |30|, |31| have been subse 



quently revised ||32|| . The fuel composition for that experiment has not been published, and 
the effect of the revision on the error bars is uncertain. Thus we do not use the Savannah 
River experiment in our fit, but concentrate instead on the more recent measurements at 
Rovno [|33|, Krasnoyarsk P3 |, and Bugey [33] where sufficient information is available. 

The thresholds {E t h) of the z/NC and z/CC reactions are 2.23 MeV, and 4.03 MeV, respec- 
tively. These relatively high thresholds, particularly for the z/CC reaction, make the yield 
more dependent on the reactor fuel composition than for the usual u e p — > e + n reaction. If 
the fuel composition is known (i.e., the fraction of fissions corresponding to 235 U, 238 U, 239 Pu, 
and 241 Pu), one can evaluate the v e flux Njj(E) in units of the number of z/ e per fission, using 



the known v e flux produced by each reactor fuel |36 



The results of Refs. p3|, are expressed as the averaged cross section (in cm /fission) 



CTfi 



N 17 {E)a{E)dE, 



(3) 



where E max is the maximum energy available in the reactor spectrum (the flux Nu(E) — > 
at E max ) . 



The results of Ref. [35| have a different normalization (in cm /u e ) 



/JT" N v (E)a(E)dE 



St" ME)dE 



(4) 



which can be converted to a fission easily (Note the erratum to that reference, however). The 
measured averaged cross sections a fission for each of Refs. |33j, |3j], [35j are listed in Table I 
in units of cm 2 /fission. 

The z/CC and z/NC cross sections were calculated to NNLO in EFT(^) and were param- 
eterized in the form 

a{E) = a{E)+b(E)L 1>A , (5) 

with the computed a(E) and b(E) tabulated in Ref. |2|. It has been noted that these cross- 
sections must be updated to include the most recent determination of qa and to include 
electromagnetic radiative corrections |J7|. To change the pion-nucleon coupling from g A = 
1.26 used in Ref. || to the most up-to-date value 1.267 ± 0.004 pBfl , we multiply a(E) 
by (1. 267/1. 26) 2 and multiply b(E) by (1.267/1.26). The inclusion of the electromagnetic 
radiative corrections increase z/NC by 1.5% and VCC by 3.5% |39| , based on the method of 
Refs. 



| 41 1 . Performing the corresponding integrals of a(E)Nu(E) and b{E)N 17 {E) and 
comparing with the experiment, we extract the parameter L\ a from each measurement and 
it is shown in the last column of Table I. The central value of L\^a varies dramatically 
between measurements. This is largely due to the fact that the term involving L^a is a 
small contribution to the cross-section, in that a change in of 1 fm 3 would produce an 
(9(1%) change in the cross-sections. The error bars listed are derived from the experimental 
uncertainties; they are dominated by the statistics of each measurement which we assume 
(for simplicity) to be gaussian in our analysis. Even though some common systematic errors 
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cr /ission(10 44 cm 2 / fission) 


T ( frr,3\ 

L 'X,A\J m ) 


f*-^ Rovno 


1.1 1 ± U.lD 


17 A _i_ i q n 
i < .4 ± lo.y 


77l\TP ^ 
WM^ Rovno 


9 71 _|_ n /1 7 
Z. ( 1 It U.41 1 


— Z.U it lo.o 


Z/CC Krasnoyarsk 


1.05 ±0.12 


-1.3 ±9.5 


l/NC Krasnoyarsk 


3.09 ± 0.30 


1.8 ±8.1 


Z/CC Bugey 


0.95 ± 0.20 


-1.5 ± 17.2 


^NC Bugey 


3.15 ±0.40 


11.1 ± 11.7 



TABLE I: Average cross-sections for charged and neutral current V e d breakup (z/CC and z/NC, re- 
spectively) and the inferred values of L\ a from each of the three reactor experiments at Rovno [33], 
Krasnoyarsk [S3], and Bugey [35]. 



are present we shall treat each of the three z/CC and z/NC as an independent determination 
of L\^a to obtain L^a = 3.6 ± 7.1 fm 3 and 3.5 ± 6.0 fm 3 , respectively. The grand average 
of all six experiments results in L\ a — 3.6 ± 4.6 fm 3 . These averages were obtained by 
summing over the individual results [L\^A)i with the error bar Si 

Zm = £%^/£-L (6) 

with the total uncertainty 

Clearly, the averaging resulted in a substantial reduction of the uncertainty. For this to be 
reasonable, the measurements must be truly independent. The total x 2 °f the averaging 
procedure is only Xtot = 2- This relatively (but not alarmingly) low value might suggest 
that there could be correlations between the z/CC and VNC results of each set of experi- 
ments. Possible sources of these correlations include: the simultaneous extraction of 17CC 
and z/NC cross-sections through the measurement of two and one neutron events for z/CC 
and, respectively, z/NC, and; the dependence on the calculated antineutrino spectra (though 
this is a small effect). It is difficult to speculate on other possible sources of correlations in 
the data. 

In addition, there are ^ 3% errors in eq. (H) due to higher order corrections. This results 
in a systematic error of ^ 3 fm 3 on L\ a- Summing up these two uncertainties in quadrature, 
we have 

Reactor = ^ ± ^ ^ (g) 

which is consistent with the constraints from tritium beta decay and helioseismology, yet 
totally independent of both. 

Given the importance of this result, one may ask whether a new high precision experiment 
of the reactor z/ e induced deuteron breakup could result in a substantially more precise 
determination of Li a- This seems unlikely, since an improvement in accuracy to ±3 fm 3 
(for example), equal to the above-stated systematic error, would require the measurement 
of the cross section afi Ssion to ~3% accuracy. Such a measurement seems rather challenging 
at the present time. 
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Finally, it should be noted that all reactor deuteron breakup data were obtained in very 
short baseline experiments. Thus, they are unaffected by neutrino oscillations, given current 
constraints []38] . 



V. CONCLUSIONS 



We have discussed how to reduce the theoretical uncertainties in the vd breakup cross- 
sections crucial for SNO's measurements of the solar neutrino flux. In effective field theory, 
the dominant uncertainties can be expressed through an isovector axial two-body current 
Li,a- We have discussed experiments that can fix L^a and have presented a constraint on this 
quantity from reactor antineutrino deuteron breakup experiments. The reactor antineutrino 
constraint alone on L^a can be translated to an uncertainty of 6-7% in the CC and NC 
reactions used by SNO, at a neutrino energy of 7 MeV. The most recent results from SNO 
verifying the existence of solar neutrino oscillations |T|], however, rely only on the CC/NC 
ratio which is insensitive to L\ t A H 0- in fact, the reactor constraint on L^a implies 
an uncertainty of only 0.7% in the CC/NC ratio at 7 MeV 0. However, a 6-7% level of 
uncertainty in cross-sections would have a significant effect on extractions of the 8 B flux. 
It is encouraging that existing analyses using the methods discussed in Section III produce 
values of Li t A which are consistent in both sign and magnitude but, clearly, more work must 
be done to constrain L 1} a to higher precision. 
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